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1.  Introduction 

co  r  -V 

.  J “ -Dege n er ale  four-wave  mixing  has  a  number  of  potential  applications  for  retro- 
reflectors,  correction  of  phase  aberrations,  laser  resonators,  signal  processing  etc.'  Most 
previous  investigations  of  degenerate  four  wave  mixing  (DFWM)  have  considered  bulk 
media  and  focussed  gaussian  beams.  In  such  experiments  it  is  impossible  to  maintain 
collimated  beams  over  long  distances  due  to  diffraction,  thus  limiting  the  quality  of  the 
DFWM  signal. 

Guided  wave  geometries  offer  some  potentially  useful  approaches  to  DFWM.2'4 
The  strong  beam  confinement  results  in  high  power  densities^and  indeed  efficient 
DFWM  has  been  demonstrated  in  both  fibers3  and  planar  waveguides.3  In  this  paper 
we  discuss  the  case  in  which  the  two  pump  beams  are  guided,  and  the  probe  beam  is 
incident  from  outside  the  waveguide  onto  the  film  surface.  Fo/'planar  waveguides. 
DFWM  has  been  discussed  previously  for  the  cases  of  all  four  waves  guided.5  and  for 
pump  beams  incident  on  a  waveguide  with  a  guided  probe  beam.6  To  the  best  of  our 
knowledge,  the  case  of  two  guided  pump  beams  and  an  externally  incident  probe  beam 
has  not  been  addressed  before.  It  is  an  interesting  geometry  for  three  reasons.  Firstly, 
the  four-wave  mixing  of  an  incident  plane  wave  occurs  effectively  over  an  interaction 
distance  of  a  micron,  decreasing  the  effects  of  distortion  etc.  Secondly.  Respite  the 
short  interaction  distance,  the  efficiencies  can  still  be  high  due  to  the  high  power  densi¬ 
ties  available  with  guided  waves,  i.e.  very  high  pump  beam  powers  are  not  necessary. 
Thirdly,  a  DFWM  signal  is  obtained  on  transmission  and  is  separated  in  angle  from  the. 
incident  probe  beam.  This  means  that  the  DFWM  signal  can  be  accessed  without  using 
a  beamsplitter  to  separate  the  probe  and  conjugate  signals. 

This  paper  is  organized  as  follows.  In  the  next  section,  we  summarize  briefly 
the  formal  aspects  of  the  calculation:  the  details  are  left  to  the  Appendix.  Section  3  is 
devoted  to  the  presentation  of  numerical  results,  and  a  discussion  of  their  significance. 


The  paper  is  summarized  in  Section  4. 
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2.  Analysis  of  the  Problem 

Hie  DFWM  geometry  analysed  here  is  shown  in  Fig  1.  The  counter-propagat¬ 
ing  pump  beams  are  both  assumed  to  be  TE„,  guided  waves,  that  is  the  electric  field 
polarization  is  the  y-direction.  (We  do  not  treat  the  TM  case  here,  although  we  do  not 
expect  anything  substantially  different  except  for  the  existence  of  a  Brewster  angle  phe¬ 
nomenon.)  The  probe  signal  is  incident  from  above  the  thin  film  waveguide  in  the 
form  of  a  plane  wave  at  an  angle  0e  relative  to  the  normal.  In  addition  to  the  usual 
reflected  DFWM  signal,  there  is  also  a  "transmitted"  DFWM  signal  whose  origin  will  be 
discussed  in  some  detail  later. 

We  now  outline  the  formal  analysis.  To  facilitate  the  understanding  of  the 
DFWM  process,  we  write  the  TE  fields  of  frequency  u  guided  by  a  film  of  thickness 
"h"  in  the  cladding,  film  and  substrate  (subscripted  c.  f  and  s  respectively)  as 

E.(r.t)  -  ±eyEs(r.u)e~iut  ♦  c.c.  -  £cnPtt\f(x)9i(ifik*  "  W°  +  c.c.  (1) 


f(x)  -  e*M.  q<.*  •  d*  -  nc*  :  x  <  0, 


(2a) 


f(x) 


Af+  eik*qf*  ♦  Af_  .**■*«.  qf’ 


n^*  -  /5*  :  h  >  x  >  0. 


Af> 


(2b) 


f(x)  -  A,e^*Mx  *  h).  q^*  -  0*  -  n,1  :  x  >  h. 

A,  -  cos(qfkoh)  ♦  ^sinlq^h),  (2c) 
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where  the  normalization  constant 


Cte  ■ 


(3a) 


+  O 


b'n  *  “  +  ifc  ♦  5K- 


(3b) 


allows  P  to  be  the  guided  wave  power  in  watts  per  meter  in  the  y -dimension.  The 
effective  index  &  is  given  by  the  usual  dispersion  relation 


t.  w  df(Qc  +  d*) 

(Qfk<>h)  "  0?"-  qfq«  ‘  (4) 

We  assume  an  external  plane  wave  incident  from  the  cladding  medium.  This 
incident  beam  undergoes  multiple  reflections  inside  the  film  and  is  partially  reflected 
into  the  cladding  (air),  and  partially  transmitted  into  the  substrate.  We  write  this  field 
as 


Einc(r.t)  -  ♦  c.c.  -  ^g^V*2  *  wt)  ♦  c.c.. 


xp  -  ncsin0c  (5) 

g(x)  -  e"*1***  ♦  Br  e'1*6***.  ice  -  nccos$c  :  x  <  0.  (6a) 

g(x)  -  B^  el*fk**  +  Bf.  e'1**15**.  -  nf*  -  xp*  :  h  >  x  >  0.  (6b) 

g(x)  -  B,  e'**1^*  *  h).  *,*  -  n,J  -  xp*  :  x  >  h.  (6c) 
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Solving  the  boundary  condition  problem  in  the  usual  way  gives 


O  „  _ _ 

"  (iCf+sX<«f+*c)exp[-2iW]  +  (ag-KfXXf-K,)  * 


- 


_ ^(Xf+ic^xp^iW] _ 

Otf+iSjXxf+KcJexpC^iW]  +  (K'-KfXiif-*,)  ’ 


B,  -  Bf*  e*w  +  %_  e*iW  . 


where  W  -  Kfk«h. 

The  calculation  of  the  nonlinear  polarization  associated  with  the  DFWM  signal 
is  relatively  straightforward.1'5  Including  the  degeneracy  factor. 


PyIL(r)  -  4e,  Q2E(x)  E.ME.ME^M . 


where  n^  is  defined  in  the  usual  way  as  a  field-dependent  refractive  index  of  the 
form  n(I)  -  n«  +  n^I  (and  I  is  the  local  intensity).  To  isolate  the  different  contributions 
to  the  DFWM  signal,  we  write 


PNl(r)  .  £  pNl(p)  «MpM  ♦  PM) 


For  the  film  region. 


Py^GQf-tf)  ■  4<oafn2Ef^TE,'^r+1®r*  • 


P^L(2qf+Kf)  -  4€0Q^n2EfCygJAf+1Bf"  ,  (10b) 

Py^Htf)  ■  8<flafn?pfCTP:A^^Af.Bf"  .  (10c) 

Py^*(+*f)  “  8<oiif a-)pfCTp:Af^Af.Bf  *  ,  ( 1  Od) 

P^*(-2qf-«f)  •  4^,^ a?pfC-rp:Af.2Bf  *  ,  (lOe) 

Pyf^'(“2q^+*f)  ■  4€oiifp?pfCTP2Af.:Bf*  .  (100 


The  dominant  wavevector  matched  terms  occur  at  p  -  hc^  for  DFWM  signal  travelling 
upwards  in  the  film,  and  at  p  -  Kf  for  the  signal  generated  in  the  downwards  direction 
in  the  film.  Multi-reflections  occurring  in  the  film  result  in  DFWM  both  on  reflection 
and  transmission  through  the  film.  Although  none  of  the  other  terms  are  wavevector 
matched,  they  also  contribute  to  the  DFWM  signals  because  the  guiding  films  are  typi¬ 
cally  a  wavelength  thick. 

A  DFWM  polarization  also  can  occur  inside  the  substrate.  Here 

pNL(-x1+2iq,)  -  4€0n,n2El(x)A*CTEJB1*eipkoh  .  (II) 

which  generates  a  wavevector  matched  DFWM  signal  travelling  upwards  towards  the 
film-substrate  interface. 

The  strength  of  the  DFWM  signal  on  reflection  and  transmission  is  calculated 
by  solving  the  polarization  driven  wave  equation  for  fields  with  wavevector  component 
-<pk0  parallel  to  the  surface,  i.e. 
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I  d* 
k^dv 


Ey(OJ.X) 


£  f  P«*(p)  eipk<,X  . 


P 


(12a) 


where 

E,(r.t)  -  ^  ey  Ey(ujc)  e^P**2  '  wt)  ♦  c.c  .  (12b) 

The  form  of  the  DFWM  signal  fields  far  from  the  film  is  given  by 

Ey(wjt)  -  x  <  0  .  (13a) 

Ey(wjc)  -  D>ek*ko(x"h)  x  >  h.  (13b) 

From  Eqn.  (6a).  and  noting  that  the  intensity  S  is  given  by  1/2  n<£e0|  Ey|  *,  the  (intensity) 
DFWM  reflection  and  transmission  coeficients  are  given  by 

i^-|Dc|*-i»fP,P.  .  (14a) 


(l  4b) 


A  variety  of  techniques  can  be  used  to  solve  Eqns.  (12)  for  R  and  T.  We  out¬ 
line  in  the  Appendix  a  technique  in  which  the  total  fields  generated  by  the  nonlinear 
polarization  source  terms  are  evaluated.  The  interested  reader  can  find  the  details 
there.  Here  we  shall  proceed  to  the  numerical  evaluation  of  R  and  T. 
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3.  Numerical  Results 


Numerical  calculations  were  performed  on  two  representative  guiO  4  wave 
film  geometries,  both  utilizing  typical  nonlinear  organic  parameters.  In  the  first  case, 
the  following  parameters  were  assumed:7  ne  ■  1.  %  -  1.6.  n,  -  1.4S.  X  -  1.06  nm. 

-  10~1‘  m2/W  and  n?F,  -  0.  In  the  second  case,  the  effect  of  a  nonlinearity  in  the  sub* 
strate  instead  of  the  film  was  examined  with  n^  -  0.  and  -  10~l*  m2/W.  Rela¬ 
tively  large  oscillations  in  the  DFWM  cross-section  were  observed  for  both  cases,  and 
two  additional  calculations  were  performed  with  a  smaller  difference  between  the  film 
and  substrate  refractive  indices,  namely  tif  -  1.6  and  n,  -  1.59. 

Figs.  2  and  3  show  rjr  and  7t.  respectively,  as  a  function  of  incidence  angle  0C. 
For  film  thickness  near  waveguide  cut-off,  the  externally  incident  beam  does  not  show 
any  resonances  across  the  film  and  the  variation  in  the  cross-sections  and  7t  is 
smooth  with  angle.  However,  as  the  film  thickness  is  increased  and  resonances  begin  to 
appear  across  the  film,  distinct  maxima  and  minima  in  the  cross-section  appear,  their 
number  increasing  with  film  thickness.  Since  the  transmitted  DFWM  signal  requires 
reflection  at  the  lower  film  boundary,  the  oscillations  are  much  larger  on  transmission 
than  on  reflection.  Therefore  it  appears  feasible  to  optimize  the  signal  for  a  given  film 
thickness  by  varying  the  incidence  angle. 

The  orijn*  of  th**  transmitted  DFWM  signal  is  fairly  obvious.  Because  reflections 
occur  at  both  film  boundaries,  the  incident  wave  forms  a  standing  wave  in  the  film. 
Hence  the  DFWM  signal  travelling  upwards  in  the  film  and  reflected  at  the  film-clad-' 
ding  boundary,  and  the  DFWM  signal  from  the  incident  wave  travelling  upwards  in  the 
film  due  to  reflection  at  the  film-substrate  boundary  both  contribute  to  DFWM  on 
transmission.  However,  because  the  reflection  coefficients  at  the  film  interfaces  are 
small,  the  DFWM  signal  on  transmission  is  much  smaller  than  that  on  reflection.  On 
the  other  hand,  it  is  "background-free"  in  the  sense  that  there  is  no  other  beam  in  this 
direction  in  space,  in  contrast  to  the  case  on  reflection  where  stray  scattering  from  a 
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beamsplitter  limits  the  signal  to  noise  ratio. 

The  variation  in  r?  with  film  thickness  at  a  45°  incident  ar;'  shown  in  Fig.  4. 
As  expected,  oscillations  again  occur  due  to  multireflections  of  both  the  incident  and 
sijpai  beams  in  the  film.  When  the  index  difference  between  the  film  and  substrate  is 
reduced,  the  amplitude  of  the  oscillations  is  also  reduced,  as  shown  in  Fig.  5. 

The  effectiveness  of  using  a  linear  film  on  a  nonlinear  substrate  is  examined  in 
Fig.  6.  Clearly  the  cross-section  decreases  rapidly  from  its  cut-off  value  and  falls 
orders  of  magnitude  below  that  obtained  on  reflection. 

We  now  compare  the  detailed  behaviour  shown  in  the  preceding  figures  with 
an  approximation  based  on  plane  wave  analysis.  It  is  well  known  that  for  plane  waves 
in  the  small  signal  limit* 


R  • 


4fLn2j  2 

X 

- 


s.s. 


(15) 


where  n  -  n«  +  naS.  L  is  the  interaction  distance.  To  a  useful  approximation  S.  - 
P./htfr  where  h#rr  is  the  effective  waveguide  thickness,  as  defined  in  Eqn.  (3b).  There¬ 
fore.  rewriting  for  the  guided  wave  case. 


R 


4#Ln 


2! 


h,frX 


P  P 


(16) 


should  provide  a  useful  approximation.  For  thick  films.  h#ff  a  h  a  L  and  the  cross-sec¬ 
tion  coefficient  becomes  independent  of  film  thickness.  Estimating  for  the  present  case, 
we  obtained  0.16xl0"lT  for  the  cross-section,  in  good  agreement  with  the  "average” 
value  of  the  oscillations  of  0.2xl0**lT  shown  in  Figs.  4  and  5.  Despite  the  fact  that 
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h  •  L  is  usually  good  to  at  best  *25%.  Eqn.  (16)  does  give  a  useful  value  for  the  cross- 
section.  In  addition  it  predicts  the  asymptotic  behav;c  -  with  increasing  film  thickness. 
Furthermore,  since  h#ff  ■*  »  at  cut-off  where  the  substrate  field  degenerates  into  a 
plane  wave  travelling  parallel  to  the  surface,  R  ■*  0  at  cut-off  can  also  be  understood. 
Finally,  it  is  clear  that  the  asymptotic  value  of  the  cross-section  with  film  thickness 
depends  primarily  on  the  nonlinearity  of  the  film,  and  not  the  details  of  the  waveguid- 
ing  structure. 

Equations  15  and  16  also  provide  some  insight  into  the  maximum  values  of  reflec¬ 
tivity  that  might  be  available.  Assuming  power  densities  approaching  damage  values, 
for  example  10  GW/cm2.  and  1  Min  thick  films,  this  corresponds  to  P  a  10T  W/m  lead¬ 
ing  to  a  maximum  reflectivity  of  0.02  (2%).  (For  1  cm  guided  wave  beam,  a  peak 
power  of  only  100  K.W  is  required  for  the  pump  beams.)  Hence,  for  this  material 
system,  pump  wave  depletion  can  be  ignored.  For  semiconductor  materials  (which  also 
exhibit  large  absorption),  reflectivities  in  excess  of  100%  should  be  possible. 

This  work  was  supported  by  the  Air  Force  Weapons  Laboratory  (F29601-87-C- 
0052),  and  the  Air  Force  Office  of  Scientific  Research  (AFOSR-87-0344). 

Appendix 

In  this  appendix  we  discuss  the  details  of  the  solution  of  the  polarization 
driven  wave  equation  (12a).  The  nonlinear  polarization  source  terms  are  given  by 
Eqns.  8-11.  The  solution  fields  can  be  separated,  for  convenience,  into  the  solutions  of 
the  inhomogeneous  wave  equation  (Py  *  0)  and  the  homogeneous  wave  equation  (Py  - 
0).  The  solutions  to  the  driven  wave  equation  do  not  by  themselves  satisfy  the  usual 
electromagnetic  boundary  conditions.  Hence  it  is  necessary  to  also  include  solutions  to 
the  sourceless  wave  equation.  For  the  cladding  and  substrate  regions  these  are  given  by 
Eqns.  (13),  and  for  the  film  by 
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Ey(w.x)  -  Dftebt'koX  +  Df.e”i*f,CoX,  x  >  h  >  0.  (A-l) 

The  structure  of  the  polarization  driven  fields  is  determined  by  the  form  of  the 
polarization  source  terms  given  in  Eqns.  (9-11).  For  the  film  region,  we  write 


Ey(r.t)  -  Yi  |Ey(p.x)  el(Htpkfl2  '  +  c.c.  .  (A-2) 

P 


Solving  gives 


Ey(PJt)  -  .  (A-3) 

Q early  there  is  a  divergence  for  the  case  p  -  fctf  which  must  be  treated  separately. 
Recalling  that  we  are  free  to  add  or  substract  solutions  to  the  homogeneous  wave  equa¬ 
tion.  we  write  for  p  -  titf 


Eytp.x)  -  5^.ipMti  -  . 


Expanding  the  exponential  within  the  square  brackets  for  small  arguments  gives 


(A-4) 


Eyfutf ,x)  -  Si  etUfk*x  . 


(A-5) 
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thus  eliminating  the  mathematical  divergence.  Furthermore,  for  the  substrate  region 
p  -  +  i2q,  and 


.  Ut.P(x-h)  (A-S) 

Similarly,  the  H^p^)  can  be  calculated  from  Maxwell's  equations. 

Equations  A- 3  to  A-6  describe  the  polarization  driven  fields,  i.e.  solutions  to 
the  inhomogeneous  wave  equation.  They  can  be  evaluated  at  the  film  surfaces,  for 
example  on  the  film  side,  to  give 


(A-7) 


pVtf1 


PyHc^iq,) 

+  4i*,q,  ♦  4q,J  ’ 


(A-8) 


(A-9) 


cttoHjih)  -  ^ 

P**r* 


PPy(P)  ipM 

pW* 


n 


+ 


+  iitfM)  ♦  PyC-Kf)*"**^-!  +  iitfkoh) 


(HC,+2iq,)Pv(Hc,) 
+  41(1,*,  +  4d,* 


(A- 10) 


Gearly  these  fields  do  not  constitute  the  complete  solution  since  the  usual  tan¬ 
gential  Ey  and  Hj  fields  are  not  continuous  across  the  film  boundaries.  Therefore,  it  is 
necessary  to  include  the  solutions  to  the  homogeneous  wave  equation,  namely  Eqns. 
(12b),  (13).  and  (A-l)  with  the  amplitudes  Dc.  Df+.  Df_  and  D,  are  adjusted  to  ensure 
continuity  of  the  boundary  conditions.  Solving  gives 


Dc-^.  (A-1 1) 

Dcn  "  «r(* ,cosW  -  iKfSinW)Ey(O)  -  d^Eyth)  +  ^cji^yh) 

-  (*fCosW  +  i^sinWjcMoI^W)  .  (A- 12) 

-  ^(Kj-MCgJcosW  -  i(*f,+icc)i1)sinW  .  (A- 13) 


and 


D 


(A- 14) 


Din  ■  -^itgEylO)  ♦  ^(KgCosW  -  iitfSinWjEylh)  -  c^H^O) 

+  c^KfCOsW  -  WgSinWjHjlh)  .  (A- 1 5) 

These  results  are  then  used  in  Eqns.  (14)  to  evaluate  the  DFWM  signals. 
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Figure  Captions 

Fig.  1  The  DFWM  geometry  analysed  here.  The  pump  beams  are  guided  by  a  film  of 
thickness  h.  A  plane  wave  beam  is  incident  from  the  cladding  side  and  DFWM  is 
obtained  both  on  reflection  and  transmission  through  the  film. 

Fig.  2  The  reflection  DFWM  cross-section  rjT  versus  angle  of  inidence  0C  for  three  dif¬ 
ferent  film  thicknesses.  Here  nc-l.  nf-1.6.  a, -1.45.  n^-lxlO"1*  m*/W  and  n2Il-0. 

Fig.  3  The  transmission  DFWM  cross-section  \  versus  angle  of  inidence  0C  for  three 
different  film  thicknesses.  Here  nc-l.  nf-1.6.  a, -1.45.  n^-lxlO*1*  mJ/W  and  n2Is-0. 

Fig.  4  The  reflection  (i?r)  and  transmission  (nt)  DFWM  cross-sections  versus  film  thick¬ 
ness  for  0c-45°.  Here  nc-l.  nf-1.6.  11,-1.45.  n^-lxlO-1*  m*/W  ad  n21l-0. 

Fig.  5  The  reflection  (ijr)  DFWM  cross-section  versus  film  thickness  for  0c-45°.  Here 
nc-I,  nf-1.6,  n,-1.59,  n2If-lxlO~‘*  mJ/W  and  n^-0. 

Fig.  6  The  reflection  (nraR)  and  transmission  (i?t*T)  DFWM  cross-sections  versus  film 
thickness  for  0c-45°.  For  the  solid  line  ne-l,  nf-1.6,  n,-1.45,  n2h-lxlO~'*  ma/W  ?«d 
n2If-0.  For  the  dotted  line  nc-l,  nf-1.6,  n,-l.59.  n2Ii-lxlO*“  m2/W  and  n2If-0. 
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